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ABSTRACT 


A potential cause of helicopter impulsive noise, commonly called 
blade slap, is the unsteady lift fluctuation on a rotor blade due to 
interaction with the vortex trailed from another blade. The relation- 
ship between vortex structure and the intensity of the acoustic signal 
is investigated. The analysis is based on a theoretical model for 
blade/vortex interaction, due originally to Widnall. Unsteady lift on 
the blades due to blade/vortex interaction is calculated using linear 
unsteady aerodynamic theory, and expressions are derived for the 
directivity, frequency spectrum, and transient signal of the radiated 
noise. The inviscid rollup model of Betz is used to calculate the 
velocity profile in the trailing vortex from the spanwise distribution 
of blade tip loading. A few cases of tip loading are investigated, and 
numerical results are presented for the unsteady lift and acoustic sig- 
nal due to blade/vortex interaction. The intensity of the acoustic 
signal is shown to be quite sensitive to changes in tip vortex structure. 
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Cartesian coordinate system, y also spanwise coordinate of 
^chapter 4 
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bound circulation on the blade (chapter 4) 
circulation in the trailing vortex (chapter 4) 
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«angle between blade and vortex, fig. 3 
^luid density 

rHifaven umber of Fourier transform in 5 , fig. 3 
-coordinate system of fig. 4 
..aerodynamic frequency (radians/second) 

-acoustic frequency (radians/second) 

.blade passage frequency (radians/second) 


INTRODUCTION 


-The sound generated by a helicopter rotor can be categorized as 
either "rotational" or "broadband" noise. Rotational noise is charac- 
^terized by a periodic acoustic signal giving rise to line spectra at the 
4)1ade passage frequency and its harmonics. Broadband noise is nonharmonic 
,.T4md has a continuous spectrum. 

-Jlotor noise is due to the combined effect of several complex acoustic 
=*unechani sms ^Sources of rotational Tioise include steady and periodic blade 
loading, and at sufficiently high blade Mach numbers-, volume displacement 
•^nd nonlinear aerodynamic effects. -Broadband noise arises from randomly 
•^varying blade forces, due for example to atmospheric turbulence. 

^nder certain flight conditions, helicopter rotor operation produces 
impulsive, highly directional noise, repeated at the blade passage fre- 
‘-squency. Commonly referred to as "blade slap," it is the predominant acous- 
tic effect when present. At least two mechanisms can be responsible for 
' ^ blade slap. One is shock formation due to local transonic flow on the 
vancing blade side. Another is the unsteady lift fluctuation on a blade 
■caused by interaction with the tip vortex trailed from another blade. This 
blade/vortex interaction is the subject of the present investigation. 

A rotor blade leaves behind a continuous vortex sheet due to spanwise 
lift variation. The inboard sheet is diffuse and does not seem to induce 
•appreciable velocities at the rotor disc, while the more intense outboard 
sheet rolls up and moves inboard, forming a concentrated tip vortex (fig. 1). 
If flight conditions are such as to cause the rotor blades to pass in close 
proximity to these tip vortices, blade slap can occur. 


- 7 - 

-Fig. 2 depicts schematically the loci of tip vortices and important 
Wade/vortex interactions for single and tandem rotors. As indicated by 
!the sketch, the essence of the blade/vortex interaction problem is that 
of a long blade passing obliquely over a vortex. For the case of small 
angles between the blade and vortex, the interaction occurs over a signi- 
-ficant portion of the blade span, so that two-dimensional aerodynamics can 
i>e used to model the problem except near the blade tip. - Comparison of 
“figs. 2a and 2b indicates that blade/vortex interactions are more likely 
rto -occur at small angles for the case of the tandem rotor, although recent 
^experimental evidence [2] indicates that ^uch interactions can be an im- 
rportant source of noise for single rotor helicopters as well. 

•**/1dnall [3] has developed a theoretical model for blade/vortex inter- 
-^ctlon. The unsteady spanwise lift distribution is computed on a two 
-dimensional airfoil passing obliquely over an infinite line vortex. This 
lift distribution is then taken as the boundary condition on a finite blade 
in the calculation of the acoustic farfield. The Widnall model has since 
•‘been refined by Chu [4] in order to incorporate sweep and convective fluid 
ef ects. 

In Widnall 's original investigation, the tip vortex is modelled as a 
potential vortex. The viscous core is taken into account by locating the 
—center of the vortex an "effective distance" below the blade, defined by 

”eff ^actual core 


Inhere r^^ is the radius of the vortex core. 

core 

The potential vortex model is convenient in that it affords an 
^essentially analytical treatment of the problem. It gives good agreement 


IS 
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^frith experiment for square tipped blades. For the case of the blade 
cutting through the vortex core, experimental results show a plateauing of 
-;both pressure fluctuation on the blade and radiated noise. This effect is 
-accurately predicted by the theory, provided has been selected 

-appropriately. 

Jn the. present study we examine the relationship between vortex 
^^tructure and the acoustic pulse, in an attempt to identify features of the 
^tip vortex that make the most damaging contribution to the radiated noise 
level. We have developed a numerical -analogue of Widnall's analysis to com- 
-|)ute the unsteady lift on the blade and the farfield acoustic signal asso- 
ciated with a vortex of arbitrary vorticity distribution. 

■We are also interested in predicting vortex structure given the dis- 
-tribution of blade loading, so that we can directly relate rotor blade 
design characteristics to acoustic performance. Recently this problem has 
received considerable attention in the study of aircraft wake turbulence, 
•-where attempts have been made to predict the trailing vortex structure in 
the wake of an arbitrarily loaded aircraft wing. Donaldson [5] has shovm 
-that the long forgotten model of Betz [6] can be successfully applied to the 
•problem. Betz assumes that the circulation and moment of vorticity of the. 
initially flat trailing vortex sheet are conserved during the rollup process 
«nd derives a relationship between the spanwise coordinate in the flat wake 
and the radial coordinate in the fully rolled up trailing vortex. 

.X)onaldson shows that complicated spanwise loading configurations gener- 
ate a discrete set of trailing vortices, such that each segment of the dis- 
tribution contained between relative maxima or minima in the loading rolls 


I 
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jjp Into a single trailing vortex. Although helicopter blade loading 
'fluctuates widely during flight, the position of maximum loading can be 
-roughly estimated at 75% of blade span^^^The loading falls off mono- 
. tonically to zero at the blade tip. Hence, we need only consider the outer 
:25% of blade span in the determination of tip vortex structure. We assume 
that the tangential blade velocity is uniform over this length. 

The Betz model is easily adapted to a numerical scheme. We apply 
••-dt here to a few different loading configurations, holding the location 
•:;and magnitude of the maximum blade loading constant for all cases. The 
ctchoice of loading configurations is somewhat arbitrary, as little is known 
-.about the precise distribution of tip loading on helicopter blades. The 
-examples have been selected to illustrate salient features of the model 
-and to suggest general guidelines for design to minimize noise. 

'Finally, we discuss some general aspects of noise generation due to 
-blade/vortex interaction and recent experimental attempts to minimize 
blade slap through tip shape modification. [2, 8] 


t 
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2. - UNSTEADY AERODYNAMICS OF BLADE/VORTEX INTERACTION 

-We wish to compute the unsteady lift distribution on a two dimen- 
sional airfoil passing obliquely over an infinite vortex filament. From 
linear aerodynamic theory, this unsteady lift can be calculated by re- 
quiring the condition of no flow through the airfoil in the presence of 
.-.the upwash field induced by the vortex. 

2.1 The Upwash Spectrum Due to Blade/Vortex Interaction 

r-The geometry of blade/vortex interaction is depicted in fig. 3. 
'Slade/vortex separation is h, angle of obliqueness A. As is standard 
-in unsteady aerodynamics, all quantities have been nondimensionalized by 
the blade velocity U and blade semi chord b„. The coordinates are fixed 
in the fluid. 

-4/e define the spatial coordinate ? measured perpendicular to the 
-cortex in the plane of the airfoil. Then given the distribution of 
circumferential velocity in the vortex, along with h and A, the upwash 
induced at the airfoil location at a given instant can be written as 
(5). Using Fourier transforms, w^(^) is represented as a distribution, 
in spatial wavenumber o, of sinusoidal gusts. 


Wo(5) 


WQ(a)e^‘^^da 


Wo(a) 


1 

2tt 


WQ(C)e‘‘"^^d5 


( 2 . 1 ) 


As the blade passes over the vortex, the upwash pattern is convected 
«long the span with nondimensional convection speed U , given by 
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u. 


1 


tan A 


( 2 . 2 ) 


This upwash can be written as w(t'), where t' is the "retarded time" 
defined by 



-We define the nondimensional aerodynamic frequency 


~ wb 

S * , the Strouhal number. 

U 


(2.3) 


(2.4) 


Here we have used the symbol in order to distinguish the aerodynamic 
frequency from the acoustic frequency that appears in Chapter 3. 

Using Fourier transforms, the unsteady upwash w(t') is represented 
•as a distribution of sinusoidal gusts of frequency S. 


w(t') 

Vf(S) 


W(S)e^^^'dS 


J_ 

2tt 


w(t' )e“^^^ dt' 


(2.5) 


We can use eqns 2.1 to express the upwash spectrum, W(S) of eqns 2.5, 
in terms of the spatial upwash distribution w^(^). An airfoil passing 
over a sinusoidal gust of spatial wavenumber experiences velocity fluc- 
tuations of frequency 

S » ccosA (2.6) 


f 


r 


•From geometry. 


cosA 


(2.7) 


•These relations are substituted into eqns. (2.5) to give tne upwash spectrum. 


»OD 

H(S) = i w^(Oe-’°«d5(^) = H^W„(o) 

« 4 » 

“ 4 


( 2 . 8 ) 


for a = 


..2-2 - Unsteady Lift Due to Blade/Vortex Interaction 


.A -gust of frequency S causes a sinusoidal lift variation along the span 


I(x,t;S) = Ljj(S,A)e 


1(k x-St) -iSt' 

* = L^(3.A)e 


(2.9) 


-■fchere Lq( 5,A) is an (as yet) undetermined complex lift amplitude. Since the 
purpose of this work -is to determine the effects of tip vortex structure 
upon the acoustic signature, we have- used the incompressible theory of 
Filotas [9] to calculate the unsteady lift due to one such oblique sinusoidal 
-gust. Compressibility effects are thereby excluded for reasons of simplicity 
.-and consistency with ref. 3. It is felt that these effects would not essen- 
tially alter the results of the present study. We are currently using the 
more complex unsteady aerodynamic theory of Amiet [15] to investigate the 
acoustic ratiation from blade/vortex interaction at higher frequencies in a 
-compressible flow. 

The Filotas theory gives the unsteady lift due to an oblique sinusoidal 


gust as 


^ 13 - 


l(x.t,S) » 

» 2iraT 

o - W(S)e1(V-5‘) (2-10) 

^where a is the Instantaneous angle of attack, 

--and T = T(ct:*A) = ^^cosA’^*''^ (complex) "lift transfer function" 

•of Pilot as, given in Appendix B. 

J^relates the magnitude and phase of -the unsteady lift to that of 
-•:the incident sinusoidal gust, analogous to the well known Sears' function 
•for the case of a parallel blade/gust interaction. T=1 at zero reduced 
■frequency (steady upwash case). As frequency tends to infinity, T 
i!>approaches zero as the infinite number of upwash oscillations become 
"anutually self-cancelling over the chord of the airfoil. 

-•Comparing eqns (2.10) and (2.9), we have for the (complex) unsteady 
lift amplitude 

LJS,A) = 2TrpU"bJ(^,A)W(S) (2.11) 


The unsteady lift at a spanwise station x is found- by integrating 
eqn (2.9) over all frequencies, as the inverse Fourier transform 


/ ® 


L(f) = 


LQ(S,A)r^^'dS 


For reference, we include the forward transform 


^(S,A) 


J_ 

2tt 


L(t')e‘*^^'df 


.S7(2.12) 




v.^ 


(2.13) 


# 
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3. -ACOUSTIC RADIATION FROM BLADE/VORTEX INTERACTION 

-Jrfe now analyze the acoustic field due to blade/vortex interaction. 
.The unsteady spanwise lift distribution previously calculated is now 
taken as a boundary condition on a blade of finite length in the solution 
of the acoustic wave equation. 

.3.1 Problem Formulation 

•4/e assume that the blade chord is smaller than a wavelength, so 
-■-=that the unsteady lift may be modelled as a line of acoustic dipoles, 
-iooving with velocity U in the y direction (fig. 4). 

The governing equation for the unsteady pressure is 


y2p 


= _L^ 

c^ 3t^ 


(3.1) 


'• Defining an effective blade/vortex interaction length s^, the boundary 
scondition may be written 


• P*(x,y,0) - p''‘(x,y,0) = L(t-x/U )5(y-Ut) 


X >s. 


(3.2) 


Dimensional variables have been introduced temporarily for the sake of 
clarity. L(t-x/U ) is the unsteady lift of eqn. (2.12), with U the 

W 

phase velocity with which the lift pattern is convected along the blade, 
from eqn. (2.2). 
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-Returning to a nondimensional formulation with reference quantities 
-Ablade velocity U and semichord b^, this boundary condition is written 


-P'(x.ytO) -p*^(x,y,Q) = L(t-x/U^)6{bQ(y-t)} Ixl^s^ 

* 0 


(3.3) 


-.where the dimensions of pressure p-force/length^ and section lift 

-J.r-'force/ length have been retained. 

* 

-3-2 Transforming the Problem 

»We -introduce the Fourier transform in x, y, and t, defining the 
--transform pair 


?p(x,y,2,t) = 


(27r)3 J 


ik x + ik V -iSt 

P(k^,k,;,S;z)e ^ dk„dk..dS 


X’ y = 


X y 


(3.4) 




p(x,y, 2 ,t)e 


-ik^x-ik^+iSt 


dxdydt 


^vhere S is nondimensional frequency, 


U 


The Fourier transform of the boundary conditions is given by 

eo S 

^ -ik x-ik v+iSt 

L(t-x/Uj5{b^(y-t)} e ^ ^ dxdydt 


AP(kx,ky,S) = 


c O' 


ji 

-« -s^ 


(3.5) 


f 
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From the properties of the delta function, the only contribution to 
^he y-integral occurs at y=t. Using 


^6(au)du = “ , 


:'--*we perform the y integration. 

• s. 




|-L(t-x/U.) e 
°o ■ 


-ik„x-ikyy+iSt 


dxdt 


— « -s. 


(3.6) 


To evaluate the t integral, we introduce the retarded time t', as 
-in eqn. (2.3). 


t' = t - 


U, 


-dt' = dt (x is held constant over the t integration) 


so that 


AP(k^,ky.S) = 


-ik X i(S-k^)(f+x/U^) 
L(f)e ^e ^ ^ dxdt' 


(3.7) 


-s. 


Regrouping, 


AP(kj^.ky,S) = 


J{C(S-ky)/U^]-k^}x' 

-s. 


i(S-kJf 

L(t')e ^ dt' 




‘ dx (3.8) 


The bracketed term in eqn. (3.8) is the Fourier transform of the 
unsteady lift in the frequency parameter S-k^. With reference to 


f 
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-eqn. (2.13), we therefore identify S-k^ as the aerodynamic frequency S, 
'Kjiving the transform of the unsteady lift as the unsteady lift spectrum 
Lq(S) (within a factor of 2tt). Since P is the Fourier transform of the 
••acoustic signal, we see that a Fourier component of unsteady lift of 
frequency S on the moving boundary will be perceived at acoustic frequency 
-S in the sound field. This is usually referred to as the Doppler shift. 
--£qn. (3.8) then becomes 


.0KP(k^,ky,S) 


2ttLq(S) 



i(S/Uc - k^)x 


dx 


(3.9) 


v-The X integration is straightforward, giving the Fourier transform of 
^he boundary condition as 


-AP(kj^,ky,S) 



sin{Sg(i/U^ - k^)} 
(S/U^ - 


(3.10) 


We now take the Fourier transform of eqn. (3.1), giving 


— - (k^^ + k 2 - k^)P = 0 (3.11) 

dz^ ^ ^ 


where the relation k=S/c has been used. 

The solution to eqn. (3.11) is written 

±zA 2+k 

■P(k^.ky,S) = A(k^,ky,S)e ^ ^ (3.12) 


where A can be complex. 
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■The argument of the exponential in eqn. (3.12) must be purely 
-complex for waves that are radiated into the farfield. This establishes 
.as-a cutoff condition 

> k ^ + k ^ for non decaying waves (3.13) 

^ Jf 


^^rom geometry k = '‘-k , so that eqn (3.12) can be written 

£ A jr 


iik,2 

-P(k^,ky,S) = A(k^.ky,S)e ^ 


(3.14) 


"^Considering only those waves that propagate outward from the boundary, 
-^0 that the quantity ±k^z is always positive, we define 

. . iA^-k ^-k z 

-for 2 > 0, P = A e ^ 


-for z < 0, 


P‘ = A'e 






(3.15) 


-=We now impose the transformed boundary condition of eqn. (3.10). 


P" - P 


2=0 


4H.L (S) 

‘’o “ 


In addition, we require that the pressure discontinuity be symmetric 
across the boundary, that is 

3P‘‘‘ _ 9p’ 

■52 ■ 'Sz , • 

2 = 0 

leading to A'*’ = -A". 


9 


(3.16) 
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✓ ■ 


vJlpplIcation of these two boundary conditions determines the constants 
A* ^and A", giving the Fourier transform of the pressure as 


P*(k^.k^.S) 



s1n{s„(i/U^-k^» 

(S/U^-k^) 


(3.17) 


taking the inverse transform tlefined by eqns. (3.4) gives the pressure 
•field as 


^^(x.y.z.t) = 


»1 

(2ir)*b 






1k x+ik v-iSt 

e ^ dk^dk dS (3.18) 

X y 


3.3 Evaluating the Farfield Pressure 

The pressure in the farfield is obtained by performing the inte- 

-grations indicated in eqn. (3,18) in the limit of /x^+y^+z^-xo. We take 

-advantage of the farfield asymptotic behavior of the integrand to carry 

--out the integrations over k and k , using the method of stationary phase 

A y 

[103. 

In general the method of stationary phase utilizes the self cancelling 
oscillation of an exponential factor in an integrand, allowing the con- 
tribution of the integrand to be neglected everywhere except in the neighbor- 
hood of one or more critical points. 

We identify the "phase" of the integrand in eqn. (3.18) as 


= z(kj^ |- + k^ ^ + /k'^-k^^-k^'^} 


( 3 . 19 ) 
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Jf the phase is a rapidly increasing or decreasing function of 

or k„, the exponential e’*^ alternates rapidly in sign with a change 
X y 

In k or k . in the limit of x,y,z-*- <», these oscillations become mutually 
X y . 

^elf cancelling, so that contributions to the integral occur only in the 

•i^lcinity of tho^e points at which the phase is rendered "stationary" with 

-respect to k and k . That is, 

^ y 

= ^ = 0 f3 20) 

3x 3y “ 10 . ^u; 


^Application of eqns. (3.20) leads to the result 




(3.21) 


-Kifhere asterisks have been used to identify and 
ky evaluated at the stationary point. 

■From geometry, k * and k * describe a wave that propagates radially 
/ ^ y 

-•^outward from a point source. Physically, the method of stationary phase 
tells us that only these waves are present in the farfield limit. 

Uhitham [11] presents the result for application of the method of 
stationary phase to multiple integrations. If we write eqn.(3.18) in the 

form 


'^p(x,y,z,t) = 


F(kj^,k^,S)e e 


-izW ik^x+ik^-iSt 


dkj^dkydS 


(3.22) 


4 J i 
»QO 


where 


21 - 


■and 


F{kj,.ky.S) 


. -1 


Vo'5) 


s1n{s^(S/U^-k^)) 

(5/Uc-k,) 


W(kx.k ,k) • -i^-kx*-ky^ 


?:::then application of Whitham's formula leads to 


4p(jc.y»z,t) = 


J F(k/.ky*.S)(|^) 

^QO 


f 

det 

a^w 

1 

ak ak. 



X y 




(3.23) 


i k *x+i k *y-i zW ( k * , k * , k ) - i J sgnW" 
e ^ y X y ^ 


-Using k * and k * from eqns. (3.21), we evaluate 
A y 


-idet 


3^W 

ak^ak^ 



k * k * 

X 


k2-k *2_b *2 


k 



and 


S - S-kj,‘ • S(1-Mf) 


where the relation k = MS has been used. 
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These expressions are substituted into eqn. (3.23) to give, after 
rsome algebra 


_p(x,y,z,t) = 


2ub, 


sin{s S[tanA(^-^^^) -M^]> 

iU{S(l-M^)} 5 n r_ 

° ^ -S[tanA(l»M^) -mJ] 


r' r" 


(3.24) 


. 1S(Mr-t) 
~ 2 MS e dS 


-«here we have substituted = 1/tanA. 

c 


^Equation (3.24) could equivalently have been expressed as an inte- 
•^^ral over aerodynamic frequency S rather than observed frequency S. The 
•latter seems the more reasonable choice because it clearly indicates that 
'•we are evaluating the pressure as an inverse transform of the farfield 
-frequency spectrum. 


p(x,y,z,t) = 


2tt 


P(S) dS 


•“We change to spherical coordinates, with reference to fig. 4. 


X = rcosSsintp 
y = rcos0costi) 
2 = rsine 


f 


S and S are now related by 


S * 5(l-Mcos6cos'l^) 


(3.25) 


:»and equation (3.24) becomes 


~p(r»6>4*,t) “ 2'jip 


peo 


.. 1 


^ Lq{S(1-Mcos6cost|j)} 


(3.26) 


-jsin{s_S[tanA(l-Mcosecos'i')-Mcosesimj^]} sine +^ 

0 

5[tanA(l-Mcos6cosii^) -Mcosesin'i^] 


:3.4 Directivity of Blade/Vortex Interaction Noise 

Using relation (3.25), we define the acoustic directivity factor 
-appearing in eqn. (3.26) as 


D(0s'i'» S) 


5in{SpSCtanA-Mcos6sin»/,^_^^^^g^^^j,^])sine 

(l-Hcosecos4,)HS(|S[taiiA-Mcosesin*/(,.„^^5gj.j,5^)] 


(3.27) 


Eqn. (3.27) is expressed as a function of aerodynamic frequency S, 

In conformity with the conventional interpretation of the directivity 
factor as a representation of the acoustic field due to a source of pre- 
scribed frequency., For a given value of S, observed frequency S varies 
throughout the field, so that writing the directivity in terms of S would 
have been misleading. 
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Using egn. (3.27), the farfield pressure of eqn. (3.26) is written 


p(r^6,tj>,t) = (1 -Mcos6cost|^)^2^ 


iL^{S{l-Mcosecosi}j)}D{e,T{),S(l-Mcos6cosi}j)} 



-1S(Mr/b„-t) 

5e ° dS 


(3.28) 


or equivalently, in terms of S, 


Ms, 


^-p(r,0,ti/,t) = ^ 


iLj 3 {S}D{ 0 ,ii;,S}S e 


i[S/(l-Mcos0cost|j)](Mr/b -t) ^ 


dS 


(3.29) 


•where r has been redefined in dimensional terms. 

The directivity factor of eqn. (3,27) can be expressed in terms of 
k, the acoustic wavenumber at the blade location. 

-Using the relation k = MS, we write 


. sin{s k[tanA/^ - cosesini|;/ 

D(e.i|;,k) = :r^- ^ 

(1-Mcos0cosi|;)2{s k[tanA/f^ - cos6sim(;/ ]} 


(3.30) 


■For the case of the source dimension much smaller than a wavelength 
(kSQ«l), eqn. (3.30) is seen to reduce to the directivity factor associ- 
ated with a convected point dipole, derived in appendix A as 
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D(e,^) » (3.31) 

(l-Mcos9cos\jj)^ 

--This expression applies in cases where the important blade/ vortex 
-interaction is confined to the tip region, as in the single rotor flow 
-^geometry sketched in fig. 1. 

-For effective blade/ vortex interaction lengths of the order of the 
•J)lade length (ksQ»l), more likely to occur for tandem rotor blade/vortex 
interactions, the directivity pattern becomes highly concentrated about 
^he combination for which 


tanA cosesln#^^,^ 

l-Mcosecosi|)^^^.^ 


(3.32) 


For these critical values of 0 and ip , the directivity factor is 
-also given by the point dipole expression, eqn. (3.31). D(e,i|)^^^^) is 
•plotted in fig. 5 for M = .53 at various values of A . (This corresponds 
to the numerical case investigated in Section 6.) The maximum value of 
^crit exceeds 50® for these cases, indicating that the acoustic 

-power is beamed strongly forward. 

3.5 The Acoustic Signal Due to a Repeated Transient 

The acoustic pressure of eqn. (3.28) is that due to a single blade/ 
-vortex interaction. For helicopter application, this transient is re- 
peated at the blade passage frequency = ma, where a is the angular 
velocity of the rotor and m the number of blades. (We note that although 
the frequency content of the transient is Doppler shifted due to motion 
of the blades, the repetition of the signal is perceived at the blade 
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passage frequency for the case of a hovering helicopter.) 
• For notational convenience, we write eqn. (3.28) as 


’P(t) = 


P(S)e 


-iSt 


dS 


^0 that 


f?(S) = 2? 


p(t)e''^* dt. 


-The repeated signal is given by 


P-(t) = I p(t-nT) 

n=-oo 


--■*»here T is the (nondimensional ) period. 


T _ 2tt jJ_ 


The spectrum of the repeated signal is 




27T 


I p(t-nT)e’^^dt = J e^’^"‘^P(S) 


n=-oB 


n=-o 


We can write 


I = ZTt I 6(ST-2imi) 


n= 


f I «(S-^) = -n^ I ) 


2tt 


II L n 

m=-« 


(3.33) 


(3.34) 


(3.35) 


(3.36) 
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so that 


Sib 


Pr(S) = -pO 


0 r 

I 5(S-fih^) P(S) 

np-oo ^ 


(3.37) 


and 


P^Ct) 


fib CO 

-XT 1 

HP* 


nb 9. f 

P(in^) e ^ ‘ 


(3.38) 


^*not,„g s^- the repeated acoustic signal is then wri 
^he Fourier series • 


written as 


P^(*’»6»«l'»t) = iO-Mcosecosip)2 ^ 


fib. 


(3.39) 


“ n.lio«Jl-Mcosecost)}D{e.^,S_^(,-Mcosecos*)}s^e^^"’'”'‘^‘’o'‘> 


f 
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4. THE BETZ VORTEX MODEL 

The Betz [6] vortex model predicts the approximate structure of 
-a 'fully rolled up lift generated trailing vortex given the spanwise lift 
nlistribution on a flat wing. ' It does not treat the intermediate stages 
-between the initially flat trailing vortex sheet and the rolled up vortex, 
-^he rollup process is assumed to be inviscid and two dimensional. 

Me use the Betz method to predict the velocity profile in the tip 
cortex corresponding to an assumed form for the steady spanwise blade 
-’loading distribution. Donaldson's results -suggest that only the vorticity 
situated between the helicopter blade tip and the position of maximum span- 
.-?wise loading rolls up into the tip vortex. We assume that the. maximum 
-loading occurs at 75% of blade span, [7] and that the loading falls off 
^jDOnotonically to zero at the blade tip, so that only the outer 25% of 
4>lade radius contributes to tip vortex structure. The tangential blade 
•velocity is taken to be uniform over this length. 

'With reference to fig. 6, the spanwise coordinate on the blade is y, 
>i»r'ch y=0 at 75% span, and y=L at the blade tip. The radial coordinate in 
"the fully rolled up trailing vortex is r. We denote the spanwise distribu- 
tion of bound circulation as rjj(y), and the radial distribution of circu- 
lation in the vortex as Fy(r}. 

The result of the Betz method is the determination of a function 
r=r(y) such that = T^{r) at corresponding values of r and y, subject 

to the constraint that r. and r, be related by three conservation relations 
for two dimensional systems: 


f 
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I. The circulation is conserved. 




drfa(y) 

dy 


dy = 


max 


5; 

dr 


dr 


(4.1) 


II. The centroid of vorticity remains at a fixed spanwise 
location. 

'VO 


7(0) - J- 


dTb(y) 

dy 


dy 


(4.2) 


III. The second moment of vorticity is conserved, 


0 = 


[y{0)-yy-^ dy 


’’max dr (r) 

r^-^ dr 


(4.3) 


"Where r is the radius within which all the vorticity 
max 

Is contained in the fully developed vortex. 

" Equations (4. 1,2, 3) are further assumed to apply piecewise, beginning 
at the wing tip, to successive portions of the sheet in toward the wing 
root: 


fyi 

dy 

L 


dy 


dry(r) 

dr 

0 


dr 


7iy,) 





d^b(y) 

dy 


dy 


L 


(4.4) 


(4.5) 
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[y. 

U(yi) 




<*^b(y) 

d> 


dy » 


dr 


dr 


(4.6) 


^ossow [12] uses the above equations to -obtain the relation between 
-r and y as 


w 

•y 

rt(y,)dy. 

4 

L 


(4.7) 


-The circumferential velocity in an axially symmetric vortex is 

41 - . . 

■!^ound from 


Ty(r) _ rjj{y(r)} 

'^e ~ Hr “ 2if 


(4.8) 


f 
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5. APPLICATION TO A NUMERICAL SCHEME 


In this section we outline the adaptation of previously derived 
-^results to numerical computation. 

5.1 The Betz Model ’ - 

Using eqns. (4.7) and (4.8) it is quite straightforward to develop 
a numerical procedure to predict the structure of the trailing vortex, 

-t^-j^iven an assumed form for the spanwise distribution of bound circulation 
fcon the blade Tjj(y). 

^Referring to section 4, eqn. {4-7) is integrated numerically from 
■y * 0 (corresponding to r = 0 in the rolled up vortex) to y = L (corre- 
responding to r„.„). For each pair of values (r, y) thus found, eqn. (4.8) 

"Is used to evaluate the circumferential velocity Vg(r). The vortex is 
.assumed to be irrotational for >'>*',j,ax* ^ potential vortex 

-'-model (Vq ~ 1/r) for this region. 

We note that Vg = 0 is always used at r = 0. (This is necessary for 
^subsequent computation.) 

' 5.2 Blade/Vortex Interactions . , 

From the circumferential velocity distribution computed as in section 
5.1 and blade/vortex separation h, we begin by computing the upwash in the 
plane of the airfoil (w^ (?) of eqns (2.1). The Fourier transform indicated 
in eqns. (2.1) is then comput6id numerically, using a fast Fourier transform 
(FFT) subroutine. The application of the FFT subroutine is described in 
-reference [13'].v 

•'Given blade/vortex interaction angle A, the upwash spectrum at the 
airfoil location W(S) is obtained from eqn. (2.8), and used in eqn. (2.11) to 
give the spectrum of the unsteady lift Lq(S). The inverse transform indi- 


f 
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cated in eqn. (2.12) is then evaluated through the use of the FFT sub- 
"Toutine to give the unsteady lift due to a single blade/vortex interaction, 
L(f). 

-Similarly for the pressure, the inverse transform of the farfield 
^spectrum is evaluated using the FFT subroutine to give the farfield signal 
-(due to one blade/vortex interaction) of eqn, (3.28). Strictly speaking, 
- 4 ue should evaluate the Fourier series of eqn. (3.29) to give the signal due 
“to a series of transients repeated at the blade passage frequency n. 
-t'itowever, numerical evaluation of eqn. (3.28) shows that the nonzero extent 
iOf the transient due to a single interaction is of much shorted duration 
.'“than the blade passage period T, so that the transient decays too rapidly 
-“to contribute acoustic energy to any subsequent pulses. Therefore, an 
^interval of the transient signal of duration T looks identical to one 
-period of the signal due to repeated interactions. 
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6 . RESULTS AND DISCUSSION 

The model Is applied to three tip loading configurations. For the 
vortex corresponding to each of these cases, numerical results are pre- 
sented for the unsteady lift experienced by a blade interacting with the 
vortex and the associated acoustic signal over a range of interaction 
geometries. General features of the noise due to blade/vortex interaction 
are considered, and results of some recent experimental attempts to 
minimize the problem through alteration of tip vortex structure are dis- 
cussed in light of the theory. 

6.1 Cases Investigated 

Inputs to the model are: 

(1) bjj, blade semichord • f 

(2) Ug^^, effective blade velocity 

(3) rjj(y), distribution of bound circulation on the blade tip 
shedding the vortex for input to the Betz model 

(4) h, blade/vortex separation 

(5) A, blade/vortex angle 

(6) Sg, effective blade/vortex interaction length 

(7) r,0 tip, coordinates of the observer, fig. 4 

(8) fi, blade passage frequency 

(9) p,c, fluid density and sound speed, taken for air at STP 
Helicopter parameters used to perform the computations were those 

of a Bell model UH-IH. Relevant specifications are listed below (from 

ref. [14]). 
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Number of Blades 
Rotor Diameter 
^Blade Semichord 
~<dvg. tip speed 
-::^vg. gross weight 
.blade passage frequency 


2 

48 ft. 
10.5 in. 
730 ft/sec 
-8000 lbs. 
•65 rad/sec 


To obtain an order of magnitude estimate of the maximum section load, 
>sinaximum loading was assumed at 15 % blade radius, with spanwise loading 
>iaken to fall off linearly from a maximum blade root 

^nd~tip. For two 24 ft blades supporting an 8000 lb helicopter, this gives 
L|jjax ® 333 Ibs/ft. Denoting maximum circulation by r^, the relation 
L * pur. from lifting line theory (with U at 75% radius) gives 
^ 238 ftVsec. 

The three spanwise distributions of bound circulation assumed for 
••the blade tip (outboard 25% span) are presented in fig. 7. is the same 
“for the three cases (elliptical, linear, and cos^ circulation distributions). 

-For the acoustic parameters, effective blade/ vortex interaction length 
Sq was estimated as 2/3 of the blade length, or s^ = 16 ft. Effective blade 
•velocity was taken to be .75U^^p, or = 585 ft/sec, corresponding to 
Mach number of .53. Calculations were performed for blade/vortex angle 
A = 0, 15, and 30®, with blade/vortex separation distance h = 0, .2, .5, 1, 
and 2 blade semi chords investigated for each value of A. 

Observer coordinates were selected to yield the most intense acoustic 
-signal for each case. (0, 4'^.^.^ corresponding to from eqns. (3.31) and 
(3.32) were computed for each value of A.) r = 1000 ft was used. 
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The cases investigated are summarized below. 

I. -Jip Load Distribution (r^ = 238 ftVsec, y = 0 at 75^ span, 

•y -« L at blade tip) 

I) elliptical r = 

II) Linear r = rQ(l - y/L) 

III) cos^ r = rQcos^(^) 

» 

JI. . Blade/Vortex Interaction and Acoustic Parameters 

I® 

Vf * 

r = 1000 ft 

•Jl* 0®, 6 = 40.7®, = 0® 

A = 15®. 6 = 38.0®, rl; = 23.3® 

A = 30®, 6 = 18®, 4' = 52.7® 

6.2 Results 

•Results are presented in figs. 8 through 27. 

fig. 8 - circumferential velocity profile in the vortices 

corresponding to each loading case 
fig. 9-17 - unsteady section lift experienced at a spanwise 
station during blade/vortex interaction 
fig. 18 - 26 - farfield acoustic signal due to blade/vortex 
Interaction 


fig. 27 


- comparison of acoustic intensities 
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6.3 Discussion 

1) Trailing Vortex (fig. 8) 

The structure of the trailing vortex, as predicted by the Betz 
theory, is profoundly influenced by the slope of the loading function at 
-the blade tip. In the case of the elliptical load distribution, for 
fwhich the loading falls off with infinite slope at the blade tip, a square- 
—root singularity in circumferential velocity is predicted at r = 0. For the 
-cos^Joad distribution, zero slope at the blade tip leads to Vq = 0 and 

d V 

0 at r = 0, and the maximum circumferential velocity occurs at a 
.^finite distance from the vortex center. The linear load distribution has 
a finite nonzero slope at the blade tip, giving a finite nonzero velocity 
at r * 0 (requiring a velocity discontinuity here). 

-Adaptation of the theory to numerical computation requires Vg = 0 
at r = 0. This acts to reduce the severity of the singularity and the 
-discontinuity encountered in the elliptical and linear loading cases. The 
-effect of viscosity in the real fluid would also moderate these abrupt 
features of the velocity profile. (Viscous effects were not considered 
•here to avoid unnecessary complication of the analysis.) 

ii) Upwash 

The distribution of circumferential velocity in the trailing vortex 
represents the upwash encountered by a blade at zero blade/vortex separa- 
tion. With an increase in separation, we see two effects: 

(1) The amplitude of the upwash signal is decreased, as circumferential 
velocity falls off with an increase in distance from the vortex 
center (except in the immediate vicinity of r = 0 for the cos^ 
loading case). 
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(2) The high frequency content of the upwash transient is diminished; 
the signal is smoothed out as the blade passes further outside 
-Abrupt fluctuations in circumferential velocity that are con- 
xentrated about the center of the vortex. 

-The blade experiences the most concentrated , sharply peaked upwash 
signal (i.e. largest high frequency content) through interaction with the 
-.trailing vortex for the case of elliptical loading. A more gradual 
variation in upwash is experienced for the cases of the linear and cos^ 
--loading distributions. The peak to peak amplitude of the upwash signal is , 
-largest for elliptical loading, less for the cos^ case, and smallest for the 
xase of linear loading. With an increase in sweep angle (A), the upwash 

C 

^signal becomes less concentrated, i.e. the high frequency content of the 
signal is diminished. 

ill) Unsteady lift (figs. 9-17) 

--The character of the unsteady lift signal can be expected to resemble 
that of the upwash signal. Accordingly, the peak to peak amplitude of the 
• unsteady lift signal is largest for the case of elliptical loading, less 
for the cos^ case, and smallest for linear loading, corresponding to the 
■ relative amplitudes of the upwash signals. An increase in sv/eep angle 
spreads out the unsteady lift signal as it does the upwash. 

The relationship between unsteady lift and upwash is governed by the 
Filotas "lift transfer function", which gives a decrease in the amplitude 
of the unsteady lift signal with an increase in gust frequency and/or 
sweep angle. We note also that the transfer function is unity for a 
zero frequency gust (steady upwash case), but introduces an aerodynamic 
•iphase shift between upwash and unsteady lift at higher gust frequencies. 


! 
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Consequently, the degree of harmonic distortion in the unsteady lift signal 
(relative to the upwash signal) is seen to .increase \^th.* an increase in the 
-Jiigh frequency content of the upwash signal. This is particularly evident 
for the case of elliptical loading at small values of h, as this upwash 
■signal is the richest in high frequency content, 
iv) Acoustic Signal (figs. 18 - 26) 

'J^e have modelled the helicopter blade as a finite line of acoustic 
-illpoles. It is not surprising, then, that analysis of the acoustic radiation 
'from a single convected point dipole affords much insight into the problan. 

-Jhe farfield sound due to a convected point dipole of amplitude 
^(point force) and sinusoidal frequency is derived in appendix A as 

“o 


where the directivity factor D(e,i|;) = ' (T-Mcosec ' os 


Two features of the dipole expression are of interest here; 

(a) The magnitude of the acoustic pressure is proportional to 
frequency, 

(b) As indicated in section 3.4, convection acts to concentrate the 
-acoustic power in the forward direction, giving rise to an in- 
crease in the value of the maximum acoustic pressure attained 
(over that of the nonconvected case). 


(a) Frequency Dependence 

The results clearly indicate the dependence of the peak acoustic 
pressure on frequency. We find, for example, that at h = 0 the amplitude 
of the upwash signal associated with elliptical loading is about twice that 
for the cos^ case, but that the acoustic signal for cos^ loading is 34 db 
-vdown. This dramatic difference (corresponding to a fifty fold decrease in 
■acoustic pressure) can be attributed to the abrupt (high frequency) 
-character of the unsteady lift signal of the elliptical case relative to the 
- -^gradual (low frequency) character of the cos^ case. Similarly, for small 
-'^values of h, cos^ loading gives a signal “that is 10 db down from the corre- 
•^ponding linear- loading case, even though the amplitude of the unsteady 
lift signal for the cos^ case is larger. 

(b) Effect of Sweep ‘ 

The introduction of sweep causes a substantial reduction in radiated 
-noise, as demonstrated in fig. 27. As already indicated, an increase in A 
-diminishes the high frequency content of the unsteady lift signal, which 
can be taken as a partial explanation for the change. 

A more subtle effect is the relationship between A and the convection 
-speed of the lift pattern through the fluid. Since the convected Mach 
number of the wave along the moving blade is M/tanA, the effective Mach 
-number of the wave through the still fluid is 

«eff = - H/s1nA (6.1) 

SO that M = sinA is the boundary between supersonic and subsonic wave 
speeds. 
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Disturbances moving at subsonic speeds with respect to the still fluid 
•radiate no sound; for M <sinA the only sound that is generated originates 
-at the ends of the blades and is due to the finite dimensions of the source. 
For this case values of (e, satisfying eqn. (3.32) do not exist, so 

that no concentration of the acoustic signal occurs. As a result, the 
radiated noise is greatly reduced. (This case was not investigated 
.•numerically. ) 

-Xonversely, disturbances moving at supersonic speeds are efficient 
-acoustic radiators, and the intensity of the peak acoustic signal increases 
•-swith Mach number. With reference to fig. 27, the diminution of the peak 
“Signal for A = 30® is due, in part, to the fact that the disturbance 
n/elocity is barely supersonic for this case (Mg^.p = 1.06). 

At large values of A it must be remembered that the two dimensional 
Tflerodynamic model used is no longer valid, due to three dimensional effects 
at the blade tip. At some point, the tip/vortex interaction becomes the 
. predominant acoustic source. 

6.- Blade-slap Minimization Through Tip Modification 

Potential methods of blade slap minimization include reduced tip speed, 
reduced disc loading, changes in blade design (area, twist, and shape) and 
alteration of blade tip shape.^^^Of these, alteration of tip shape is of 
■particular relevance to the present study, as this technique entails changes 
in tip vortex structure that minimize noise. 

The most significant result of this investigation with regard to tip 
design is that the slope of the spanwise loading distribution at the blade 
tip strongly influences tip vortex structure and the associated intensity of 
noise due to blade/vortex interaction. This result is consistent with recent 


experimental findings [8] that show substantial reduction in blade slap 
Intensity through the use of a tapered ("ogee") blade tip. (fig. 28a) 
Mother technique that has yielded positive results [2] makes use 
of a "sub-wing" tip (fig. 28b). This design apparently divides the tip 
•vortex into two less intense twin vortices that destructively interact 
to produce a diffuse trailing vortex. However, more quantitative acoustic 
data is needed to judge the efficacy of this method. 
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7 . CONCLUSIONS 

'^lade/ vortex interaction is an important source of helicopter blade 
slap. Several experimental studies indicate effective blade slap mini- 
mization through modification of tip design in an effort to alter vortex 
structure. 

-Theoretical analysis of blade/vortex interaction, through application 
of linear unsteady aerodynamics and the acoustics of a moving source, is 
^ble to account for the essential role played by many of the parameters 
■iithat govern the occurrence of blade slap. 

~ Jhe -numeri cal model developed to extend the analysis to tip vortices 
:of -arbitrary structure constitutes a powerful investigative tool. The 
-jnethod affords useful insight into the relationship between the character 
of the tip vortex and the associated noise due to blade/ vortex interaction. 

The intensity of the noise was found to be quite sensitive to changes in 
•tip vortex structure. 

This model is applicable to blade/vortex interactions at small 
ob’ique angles, for which the interaction occurs over a large portion of 
■the blade span. Such interactions are probable for tandem rotor helicopters. 
For single rotors, blade/vortex interactions are likely to occur at large 
angles, for which the interaction is confined to the blade tip region. 
Additional theoretical modeling is required to account for the three 
dimensional aerodynamic effects associated with this case. 

Other effects the model is not equipped to handle include blade or 
self-induced motion of the vortex filaments during blade/vortex interaction, 
-curvature of the vortex filaments, and blade rotation. These problems deserve 
future study. 


The aerodynamic theory used to model ijlade/vortex interaction 
•strictly applies only to the case of a low frequency gust in an in- 
compressible fluid. Making the appropriate low frequency assumption that 
the blade chord is smaller than a wavelength, we neglect chordwise pressure 
• variation in the unsteady lift calculation and model the blade as a line 
s4)f dipoles in the analysis of the acoustic field. At higher gust fre- 
-.quencies or higher Mach numbers, compressibility effects become important, 
««nd a more powerful aerodynamic theory such as that due to Amiet [15] is 
,-Tequired. The simpler aerodynamic model used here is sufficient to account 
•for the essential details of the relationship between vortex structure and 
^he acoustic pulse. The additional refinement obtainable from consideration 
.of compressibility and high frequency effects was felt to be unnecessary 
for this initial study. The extension to include the full effects of com- 
-pressibility is currently underway. 


—APPENDIX A: Sound Field Due to a Convected Point Dipole 


The calculation of the sound field due to a convected point dipole 
follows the analysis of section 3 (for a moving line of phased dipoles). 
The simpler boundary condition, however, gives the farfield preslure sig- 
-nal in integrable form for the case of harmonic oscillation. 

Jrfe consider a point force of magnitude f(t), oriented normal to the 
plane, moving v/ith velocity U in the y direction. 

The governing equation for the unsteady pressure is 

\ 

V2p = i_ iffi. (Ai) 

c* 3t^ 

•^ith boundary conditions 

Ap = p(x,y,0‘) - p(x,y,0'^) = 5(x)6(y-Ut)f(t) (A2) 


Wt introduce the Fourier transform in x,y, and t, defining the transform 

pair 


p(x.y,z.t) = 


1 


(27t) 


P(k^.ky.a);z) = 


, i(k x+k v-ut) 

P(k^.k^.w;z)e ^ ^ dk^^dk^dcj 


-1(k^x +k^-o)t) 

p(x,y,z,t)e ^ ^ dkdydt 


The Fourier transform of the boundary condition is given by 

-1(kj^x+k^-o)t) 


AP(k^,k^,o)) = 




. J 

’ « 


6(x)6(y-Ut)f(t)e 


dxdydt 


(A3) 


(A4) 
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From the properties of the delta function, the only contribution to the 
y Integral occurs at y = Ut, giving 




-Ik x+i (ai-k U)t 
f(t)6(x)e ^ ^ dxdt 


We define the frequency w = co-k^U, recognizing that a frequency w 
located on the moving boundary is perceived at Doppler shifted frequency 
0 ) "in the field. Eqn. (A5) is then written 


-ik X ffoo iwt ' 

'AP(k^,k„,w) 5(x)e ^ J f(t)e dtl dx 


Defining the transform pair 


-f(t) = F(S)e"*“*dS 


FG) “ ^ f(t)e’“‘dt 


'we Identify the bracketed term in eqn. (A6) to be 2 ttFQ), 


-ik X 

AP(k ,k , 0 )) = 2ttFQ) 5(x)e ^ dx 

X Jr 


Contribution to the x integral occurs at x = 0, giving the Fourier 
transform of the boundary condition as 


1 
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AP(kx*ky.w) = 2TrF(w) 


(A9) 


where w = <»i-kyU 


Taking the Fourier transform of eqn. (Al) leads, as in section 3, to 
the expression 


^ ^ ‘-k ‘ 

P"(x.yt2.t) = A"(k^,k„.(i>)e ^ 


x’y 


(AlO) 


: * ; 

I: 


.^Applying the boundary condition of eqn. (A9) and imposing the 
•additional requirement 


- 

32 


3P* 

3z 


2 = 0 


(All) 


-'gives the constants A~ = + irF{(ij). 

The unsteady pressure is then given by the Fourier integral indicated 
xtn eqn. (A3), 


p(x,y,z,t) = - -~ 




dk^.dk do) (A12) 
X y 


i ; 


As in section 3, the k and k integrations in the asymptotic farfield 

X y 

-limit can be evaluated by the method of stationary phase, giving the 
farfield acoustic pressure as 



p(x,y,z,t) 
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IFM 

r* c 


-1u(t-r/c) 
e du 


Xhanging to polar coordinates, with reference to fig. 4 , 
X * rcos0sint|» 
y = rcos0cosi|i 
X = rsin0 
gives 


«T)(r,0,Tj>,t) 


•9«here 


1 

4iTr 


oo 

iF(o))sin0 ^ 


J.OD 

w = (ij(l-Mcos8cos4i) 


-1u(t-r/c) 
e d(i> 


(A13) 


(A14) 


For the case of a simple harmonic dipole of frequency 


X(t) = F^e 




(A15) 


"the transform is given by 


F(u)) = 


1 

2tt 




!o 

2ir 


® i(a)-(o.)t 
e ° dt 


FoS(u-a>o) 


(A16) 


Using (o = “/(l-Mcos6cosij/),eqn. (A14) then becomes 


r » 


p(r.0,tj;.t) = ^ 


6(a)-w ) 


s1n0 w p”^^l-Mcos6cosip^^~*^^^^^~ 


° (1-Mcos0cosiji) 


2 c 


(A17) 
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Contribution to the integral occurs at w ® giving the unsteady 
-pressure as 


p(r,0.t;;,t) 


Is. !!o O(e,<,)e"’*’-M«sscosif 

2Trr c 


-where the directivity factor is defined by 




\ 


(A18) 


(A19) 
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■APPENDIX B: Lift Transfer Function of Filotas Theory 

The lift transfer function is used to calculate the unsteady lift 
-due to an oblique sinusoidal gust impinging on a two dimensional airfoil 
-in Incompressible flow (fig. Bl). For a gust of nondimensional v/ave- 
^-inumber a = 

A)(l-«- I sinA) 

-exp{-ia[cosA- ^ ]} 

l+2ir0(l + Y sinA) 


T’Ca.A) 


1 +Tra(l+cos^A.+ nosinA ) 


APPENDIX C: Computer Program 


ORIGINAL PAGE IS 
OF POOR QUALITY 



COMPLEXES CU048) «W(204d) ,CL^T«W0>L0*P 
COMPLEX+8 CMPLX,CEXP 

REAL L,LAMBDA,HACH,VTHET (1024) ,VR (1024) 

DIMENSION GBAF1 (512,7) 

DIMENSION UU(514) ,V1 (514) ,V2 (514) , V3 (51 4) *V4 (5^4) , V5 (SU) , V6 (514) 
COMMON W,C,VR,R, THETA, PSI, LAMBDA, OHEG A, GAMK A,G AMMAO, M ACH,RHO,SO, 

6 UU,U,H,P1,UFL0, 

6 COSLAM,TANLA«,COTHET, COSPSI, 

& DS, DEBAR, DT, DTP, DELZ, DR, 

6 N,N2,N4 

C VARIABLES USED IN MAIN PROGRAM: 

C BLDS; NO. OF ROTOR BLADES 

C BO: BLADE SEMICHO RD (FE RT) 

C CL:UNSrEAD]f SECTIOtN LIFT COEFFICIENT PROM PILOTAS TMBORY 

C D1REC:DIRECTIVITY FACTOR, EON. 3.27 

C GAMMAO:CIRCULATIO.N (FT**2/SEC) 

C GAHhA:NONDIHBNSIONAL CIRCULATION (GAMHA=0AMMA0/B0*U) 

C GC:GRAVITATIONAL CO NST ANT (LBM*FT/LBF*SEC**2) 

C H:VORTEX TO BLADE DISTANCE IN SEHICHORDS* SEE FIG 3 • 

C LAMBDA: ANGLE BETWEEN BLADE AND VORTEX, FIG. 3 (RAD) 

C MACH: BLADE MACH NUMBER 

C OMEGA: ANGULAR FREQUENCY OF ROTOR (RAD/SEC) 

C (CORRESPONDS TO ALPHA OF TEXT) 

P:PKE3SURE (LBF/FT*-*2) 

C ■ R, THETA, PS1:C00RDINATE SYSTEM OF FIG. 3 
.C SO: BLADE INI^RACTIpN -LENGTH IN SEMICHORDS 

C • U:BLADE VELOCITY (-^T/SEC) ' 

C R1I0:FLUID DENSITY (LBM/|FT**3) 

C SBAH: STROUHAL FREQUENCY (AERODYNAMIC FREQUENCY) 

C S: DOPPLER SHIFTED STrOUHAL FREQUENCY (ACOUSTIC FREQUENCY) 

C UFLO: USED IN UNDERFLOW CHECK' 

PI=3. 1415926 • ■ 

B0=10.5/12. 

S0=16./B0 

R=1000. 

DLGRAD=2.0*PI/360.0 • 



LAMDDA»1&. 0*D£GRAD 
THETAs=38. ♦DEGRAD 
PSI=23.3*DEGRAD 
U=505. 

KACH=U/1100. 

OMEGA=32.5 

GAHMA0=230, 

GAMMA=GAMMAO/ (B0*0) ’ . ' 

GC=32.17 
RUO=. 077/GC 
COTHBT=COS (THETA) 

COSPSI=COS (PSl) 

COSLAK=COS (LAMBDA) 

TANLAM=TAN (LAMBDA) ‘ ; 

K=2048 

N-2=H/2 

N4=H/4 

DO 100 M=1,N 
C (K)=CMPLX (0.0,0,0) 

100 COHTINUE 
TEN=10.0 

UFLO^TEN** (-40.0) 

Y0=0.0 

C H=0.0 

g ♦ TO AVOID CONFUSION DOE TO THE NUMBER OF DIFFERENT CASES * 

C * INVESTIGATED, SUBROUTINE CALLS AND PLOTTING ROUTINES HAVE NOT * 

C ♦ BEEN INCLUDED HERE. THE FOLLOWING SEQUENCE WOULD RETURN THE ♦ 

C ♦ ACOUSTIC SIGNAL IN COMPLEX ARRAY C: ♦ 

C ♦ CALL BETZ(NLOAD) (SEE SUBROUTINE BETZ) ♦ 

C ♦ CALL TRNSFM ♦ 

C * CALL SOUND ♦ 

C SUBSTITUTION OF 'CALL LIFT* FOR 'CALL SOUND* WOULD RETURN THE ♦ 

C ♦ UNSTEADY LIFT IN COMPLEX ARRAY C ♦ 

Q ^ ^ if *¥***** 4*^* ^****4l ******* 1ft M********tf1f 


nonooonoooonooonoooooo 


SUBROUTINE BETZ(NLOAD) . , 

BETZ COMPUTES THE ClHCUMFERENTlJlt VEIoCITY DISTRIBUTIOM IN THE 
TRAILING VORTEX USING THE BETZ ROLLUP MODEL OP CHAPTER 4, AND 
EVALUATES THE VERTICAL VELOCITY H AS A FUNCTION OF Z, THE 
SPATIAL COORDINATE PERPENDICULAR TO THE VORTEX. 

(Z CORRESPONDS TO XCI IN TEXT, FIG 3) 

SAMPLING INCREMENT. IS DELZ, OVER INTERVAL VRTSEP 

H(l) CORRESPONDS TO Z=0. W(t,...1024) CONTAIN VALUES FOR POSITIVE 
Z, W (1025, .. .2048) VALUES FOR NEGATIVE Z. VALUES ARE STORED IN 
ORDER OF INCREASING Z 

COMPUTATION PROCEDURE: 

EQN 4.7 IS INTEGRATED NUMERICALLY AND USED WITH EQN 4.8 TO COMPUTE 
A SET OF N/2 CORRESPONDING VALUES RD (M) AND VCOH (M) , THE RADIAL ' 


COORDINATE AND CIRCUMFERENTIAL VELOCITY RESPECTIVELY IN THE 
ROTATIONAL PORTION OF THE TAILING VORTEX. VCOR IS THEN SAMPLED 

AT INCREMENT DR OVER N/2 POINTS TO COMPUTE VR, THE w 

CIRCUMFERENTIAL VELOCITY PROFILE IN THE TRAILING VORTEX. OUTSIDE • 

THE ROTATIONAL REGION (R>RD(N/2)) VR IS PROPORTIONAL TO 1/R. 

UPWASH H IS COMPUTED FROM GEOMETRY, 


NLOAD is an INDEX IDENTIFYING THE TIP LOADING CASE INVESTIGATED. 
NLOAD=1 FOR ELLIPTICAL LOADING, 2 FOR LINEAR, AND 3 FOR COS**2. 
COMPLEXES C ( 2040 ) , W ( 2048 ) ,CL,T, W0,L0,P 
• COMPLEXES CMPLX-,CEXP 

REAL L, LAMBDA, MACH, R5(1024), VCOB(1024), ’ ' 

6VTHET.(1024) ,Vtf (.1024) ' 

INTEGER NLOAD • . • • • ' ’ r 

COMMON W,C,VRVR,TUETA',PSI, LAMBDA, OHEGA,GAMMA,GAMHAO,MACH^RHO, SO, 
G BO,u,H,Pi,uFLO, ! ' ’ 

G COSLAM, TANLAM,COTHET,COSPSI, ' ' ; 

G DS,DSDAR,DT,DIP,DELZ, DR, • ^ 

G N,N2,N4 • ■ 

BLDS=2. . ■ l.-i-.. •' 

VRTSEP= (2.0#PI/(BLDS^OMEGA^B0)) ♦U*COSLAM‘ -■ V . 
VBTSEP=VRTSEP*4. ‘ 



o o . o o o o 


■ 

i * 

i 

1 . 

I . . 

i 

DELZ=VRTSEP/PLOAT (M) 

• DR=D£LZ/1.95 

L=15./B0 
D5f= (L/2.0) /M2 
SUM=0. 0 
RD(1)=0.0 

VCOR(1)=0.0 : 

YY IS THE SPANHISE COORDINATE. (CORRESPONDS TO Y OP TEXT) TY*0 AT 
15% SPAN, YY=L/2 AT BLADE TIP. (IN TEXT, Y*L AT BLADE TIP) 

VCOa(H) IS THE CIRCUMFERENTIAL VELOCITY COPHESPONDING TO 
RADIUS RD(H) IN THE TRAILING VORTEX. (VTIIETA, H USED IN TEXT) 

DO 100 H=2,N2 
FM=«“ 1 

YY= (L/2. 0) - (FM*DY) 

GO TO (14, 15,16) ,NLOAD 

14 GG= (YY/ (L/2. 0) ) **2 
GAM=GAMMA* ( (1 . 0-GG) ♦*0.5) 

GO TO 17 

15 GAH=GA!1HA* (1.0- (YY* (2.0/L) )) 

GO TO 17 

16 GG= (YY + PI) /(L/2.0) 

GAM=0.'5*GAMMA* (1 .O+COS (GG) ) 

17 CONTINUE 
SUK = SUM+ (GAM*DY) 

RB (M) =ADS (SUM/GAM) 

VCOR(M)=GAM/(2.0*PI*RB(M)) • 

100 CONTINUE 

CORAD IS THE RADIUS HITUIN NHICH ALL VORTICITT IS . 

CONTAINED. (RMAX OF TEXT) •' • 

CORAD^RB (N2) 

VR(1)=0.0 ■ • ■ , 

DO 400 M=2, N2 . -.1. • 

FK = M-1 ■ :■ 

RP=FH*DR ... • ■ 

IF (RP.GT. CORAD) GO TO ’700 ’ \ ‘ 

DO 500 J=2,N2 ... 


noon 


IP (RD(J).GT.RP) 60 TO 
IP (J.EQ.N2) GO TO 600 
500 CONTINUE 
600 VR (M)=VCOH (J) 

400 CONTINUE 
700 CONTINUE 

DO 800 K=M,N2 

FK=K-1 

RP=FK*Dn 

VR(K)=+uAHHA /(2,0*PI*RP) 

800 CONTINUE 

DO 250 1=1, N 
W (I)=CilPLX (0. 0,0.0) 

250 CONTINUE 

DO 900 H=2,N2 

FH=M-1 

Z=F«*DEL2 

RZ= (Z**2 + n**2) ♦♦O.S 

POS=RZ/DR 

K=IFIX (POS+0. 5) 

IF (K.LT.2) K=2 
IF (K.GT.N/2) K=N/2 
ETA=ATAN2 (H,Z) 

C VRIVEL IS VERTICAL VELOCITY 

VRTVEL= 4-VR (K) *COS (ETA) 

IF (K.EQ.N2) VRTVEL=+GAMMA /<2.0*PI*Z) 

TAILI.AND TAIL2 ARE THE POSITIVE AND NEGATIVk "TAILS” OP THE 
UPWASH PROFILE THAT APPROACH ZERO AS Z APPROACHES INFINITY. THEY 
ARE USED IN THE "ALIASING" STEP (DESCRIBED IN REF. 13), NECESSARY 
FOR SUBSEQUENT COMPUTATION OF FAST FOURIER TRANSFORMS 
TAIL1=0.0 
TAIL2=0. 0 
DO 150 1=1,41 
FJ=I-21 

IP (FJ.EQ.O) GO TO 150 
Z1= <FJ»VRTSEP) +Z 


1 

f 

I 


22*(PJ*VRTSEP) -Z 

TAILl=TAIL1+GAIiaA/(2.0*PI^Zl) 

TAIL2=TAIL2*GAMHA/(2.0*PI*Z2) 

150 CONTINUE 

W (M) = CMPLX (-VHTVEL^0.0)*CMPLX(-TAIL1,0.0) 

H (N-M+2)=CMPLX (+VRTVEL,0. 0) +CMPI.X (-TAIL2^0«0) 
900 CONTINUE 
RETURN 
END 



n o n n o o n 


SOBROOTIMB tiPT ■ ( ■ i 

SOBROUTINii LIFT COMPUTE^ THE UNSTEAbY tlPT OW THE blADE BY 
EVALUATING THE INVERSE FOURIER TRANSFORM OP LO, THE COMPLEX LIFT 
AMPLITUDE. LO IS EXPRESSED AS A FUNCTION OF AEEODVNAMIC FREQUENCY 
SOAR. LIFT IS EXPRESSED AS A FUNCTION OF RETARDED TIME TP, 

WHERE TP=T-X/UC. UC IS THE PHASE VELOCITY OP THE LIFT • 

PATTERN, UC=U/TAN (LAMBDA) 

VALUES ARE RETURNED IN COMPLEX ARRAY C. 

COMPLEX + 8 C(2048) (2048) ,CL,T, W0,L0, P,V ' 

COMPLEXES CHPLX,CEXP 

REAL L, LAMBDA, «ACH,VTHET (102 4) ,Vfl (1024) 

COMMON K,C,VR, R, THETA, PSI , LAMBDA, OMEGA, GANNA, GAHNAO,HACH,RUO,SO, 

C bO,U,H,PI,UFLO, 

0 COSLAM, TANLAH,COTHET,COSPSI, 

6 DS,DSBAR, DT, DIP,DELZ,DR, 

& N,N2,N4 


177 

DO 177 K=1,N 
C (M)=CMPLX (0.0, 0.0) 
CONTINUE 



c 

A=2.0*PI*RHO* (U^*2) *80 
DTP=(2.0*PI)/(FLOAT(N) *DSBAR) 
TP0=DTP*FLOAT (N) /2 . 

DO 100 M=2,N2 
FM=M-1 

DO 200 1,2 


" ■ ■ ' 


IF (J.EQ.2) FM=-FM 

i • ; 



SBAR=FM*DSBAR 

/ ’ ' * 


C,- H , 

CL=T(SDAR/C0SLAM,LAMBDA) 



r' 

IF (J.EQ.l) W0=W(M)'v‘. 

< ' 


JO 

IF (J.EQ.2) H0 = W(N-M + 2) 


“•V * 

> 

LO=A*CL*WO 

\ * 

... 

r^, 

H-’; ' 

v< Uj 

LO=LO*CBXP(CMPLX (0.0,SBAR*TPO)) 
IF (CABS(LO) .LT.UFLO) GO TO 500 


■* ■' 

\ 


IF (J.EQ.l) C(M)*L0 . 

IF (J.EQ.2) C(N-M+2)=L0; 
200 CONTINUE 


100 CONTINUE 

C FOUNT IS A FAST POURIEti TRANSPOHM SUBROUTINE 
500 CALL FOURT(C,N,-l 1,-1, + 1, WORK) 

DO 300 M=1,N 
C (M)=DSBAR*C(M) 

300 CONTINUE 
RETURN 
END 


n rvo n n n o 


I 


SOBROOTINB SOUND , ' . . ‘ . 

SUBHOOTINE SOUND COMPUTES THE FARFIELD PRESSURE StCMAt BY 
EVALUATING THE INVERSE FOURIER TRANSFORM OF P* THE PHEQUEHCY 
SPECTRUM OP THE FARFIELD SIGNAL. P IS EXPRESSED AS A FUNCTION OP 
OBSERVED FREQUENCY 3. THE FARFIELD PRESSURE IS EXPRESSED AS A 
FUNCTION OF TINE,T. OBSERVED FREQUENCY S IS RELATED TO 
AERODYNAHIC FREQUENCY SBAR BY S3AR=S* ( 1 -M ACH^COS (THETA) ^COS (PSI) ) 
VALUES ARE RETURNED IN COMPLEX ARRAY C 
COMPLEXES C (2040) ,W (2048) , CL, T,RC,L0,P,V,VANLYT 
COMPLEX^O CMPLX,CEXP 

REAL L, LAMBDA, MACH, VTUET( 102 4) ,VH (1024) 

COMMON W,C, VR, R, THETA , PSl , LAMBDA , OMEGA, GAMM A« GAMMA0« HACH#RiiO, SO, 

6 UO,U,H,PI,UFLO, 

& COSLAM, TANLAN,COTHET, COSPSI, 

6 DS,DSBAH,DT,DTP,DELZ,DR, 

& N,N2,N4 

DO 177 M=1,N 
C (M)=CMPLX (0.0,0.0) 

177 CONTINUE 

A=2.0*PI<^RHO* (U**2) *80 
B= (MACH*S0) /(2.0*PI*R) 

X=MACn*COTHET*COSPSI 

Y=1.0-X 

DS=DSBAR/Y 

■ DT=(2.0*PI)/(FLOAT(N)*DS) 

T0=DT*FLOAT(N)/4. • > • 

DO 400 M=2,N2 . ' 

FM=M-1 . •• ■ 

DO 500 J=1,2 . ■ . ' ■ ■ 'V ■ 

IF (J.EQ. 2)- FM--FM 
s=FM»Ds • 

SDAR = FM*DSBAR , 'i- 

CL=T(SOAR/COSLAH, LAMBDA) - ' 

IF (J.EQ.1) W0=H(H) 

IF (J.EQ.2) H0=H(N-M<*2) • 

L0=A»CL*H0 #r'.V 


E=SO*S* ( (TANLAM*Y) -X) 

DIREC=SXN (THETA) *SIN (E)/fi 
P=D*LO*S»DIREC*CMPLX (0.0, 1*0j 
P=P*CEXP(CHPLX (0.0,S*T0) ) 

IF (CADS (P) .LT.UFLO) GO TO 700 
IF (J.EQ.1) C(M)=P 
IF (J.EQ.2) C(N-M + 2)=P 
500 CONTINUE 
400 CONTINUE 

C FOURT IS A FAST FOURIER TRANSFORM SOBROOTXMB 

700 CALL FOURT (C,N,+ 1,-1, <•!, WORK) 

DO 600 tt=1,N 

C(M)=DS<'C(«) * 

600 CONTINUE 
RETURN 
END 


cn o C*» c» 


c 

c 

c 

c 


c 


100 


200 


SOBROOTINB TRNSPM , 

TRNSPrt COMPUTES THE POORIBA ■tfiAN^PORM OP VERTIcAL VELOCITY 
H IN TERMS OP WAVENUMBER SlUMA, THEN EXPRESSES THE TRANSPORH 
AS A FUNCTION OP AERODYNAMIC PREQUENCY SBAR, 

WHERE SDAR=SIGMA*COS(LAMDDA) 

COMPLEXES C (2048) ,W (2048) ,CL,T,H0,L0,P 
COHPLEX*8 CMPLX,CEXP 

REAL L, LAMBDA, MACH, VTI1ET(1024) ,VR (1024) 

COMMON W,C,VR,R,TH£TA,PSI,LAMBDA,OMEGA,GAMHA,GAMHAO,HACU«RHO»SO 
B0,U,U,PI,UFLO, 

C0SLAM,TANLAM,C0THET,C0SPSI, 

DS,DSBAR,DT,DTP,DELZ, DR, 

N,N2,N4 

FOURT IS A FAST FOURIER TRANSFCHM SUBROUTINE 
CALL FOURT (H, N,f 1,-1, + 1, WORK) , 

DO 100 H=1,N 

W (H) =B (M) ♦ (DELZ/ (2. 0*PI) ) 

CONTINUE 

DSIG= (2,0»PI) /(FLOAT(N)*DELZ) 

DSBAR=DSIG*COSLAM 
DO 200 M=1,N 
W (M)=W (M) /COSLAM 
CONTINUE 
RETURN • 


o o o on 


PONCTIOM P«IN(C,tl) 

FHIN IS THE HlMIMOn BEAL VALUE CONTAINED IN A COMPLEX ABRAX C OP 

LENGTH N , 

COKPL£X*8 C(N) 

P=0.0 

DO 100 M=1,N 
Q=flEAL(C(M)) 

IF (Q.LT.P) PsQ 
100 CONTINUE 
FMIli=P 
PETURN 
END 

FUNCTION FHAX(C,N) ’ ' 

FHAX IS THE HAXinUH BEAL VALUE CONTAINED IN A COMPLEX ARRAY C OP ’ . 

LENGTH N ' 

COHPLEX*8 C(N) S 

P=0.0 • 

DO 100 M=1,N 
Q=REAL(C(M)) 

IF (Q.GT.P) P=Q 
100 CONTINUE 
FKAX=P 
RETURN 
■ END 


C> O c: CO 


SUDROUriNE DB 

DO COHPUi'ES THE INTENSITY OF THE ACOOSTIC SidNAL IN DECIBELS 
COMPLEX*a C(20«8) ,H(20ll8),CL,T, H0,L0 
COMPLEX+8 CMPLX^CEXP 

REAL L^LAMDDA, MACH, VTHET( 1024) ,VR (1024) 

COMMON W,C,VR,R, THETA, PSlfLAHDDA, OMEGA, GAMMA, 6AMHA0, MACH, RUO, SO 
B0,U, H, PI,UFLO, 

COSLAM,TANLAH,COTH£T,COSPSI, 

DS, DEBAR, DT, DTP, DELZ, DR, 

N,N2,N4 
DELP=FMAX (C,H) -FHIN (C,N) 

P= (DELP/2.0) *47.88 

PR£F=0. 00002 

DBL=20. 0*ALOG10(P/PHEF) 

WRITE (6,99) H,DBL 
FORMAT (F10,2,F10.2) 

RETURN 

END 


COMPLEX PaNCTIOM T (SIGmX, LAMO bA) 

T IS THE LIFT TRANSFER FUNCTION GIVEN BY PILOTAS THEORY 
SIGMA IS A DIMENSIONLESS HAVE NUMBER, SIGMA»SBAR/COS (LAHBOA) 
REAL LAMBDA 
COMPLEX+8 £ 

COMPLEXES CNPIX,CEXP 
S1NLAM=SIN (LAMBDA) 

COSLAH=COS (LAMBDA) 

PI=3. 1415926 
SIGABS=ABS (SIGMA) 

A=1.0+ (0.5*SINLAM) 

B=PI* ( (PI/2. 0) -LAMBDA) ♦A 
C=1.0+( (2.0*PI«SIGABS)*A) 

PiJI=COSLAM- (B/C) 

D=1.0* (C0SLAM**2) ♦ (PI*SIGABS*SINLAM) 

F=1 .0+ (PI*SIGABS*D) 

G=F**0.5 

E=CEXP(CMPLX(0.0,-SIGMA+PHI) ) 

T=E/CMPLX (G,0.0) 

RETURN 

END 


Figure 1. Generation of Helicopter Rotor Tip Vortices (after ref. 4) 
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Figure 3. Geometry of Blade/Vortex Interaction (after ref. 3) 
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Figure 4. Geometry of Acoustic Model 









Figure 6. Geometry applied to the Betz vortex rollup model. 

Maximum loading is assumed at 75% span, so that 
vorticity shed from the outboard 25% span. rolls up 
into the trailing vortex. 
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if-lgure 28a. Comparison of farfield acoustic signatures for ogee(fig.) 

jmd square-tipped rotor blades. Tests were conducted v/ith 
a Hughes model UH-IH helicopter in level flight at a 
‘-^nominal 61 meter altitude. Microphone is located at ground • 
•.level in line with the flight path. (ref. 8) 



Figure 28b. Sub-wing Tip (ref. 2) 
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-Figure Bl. Infinite Span Ai 
(after ref. 4) 
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